Modeling temporal networks using random itineraries 
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We propose a procedure to generate dynamical networks with bursty, possibly repetitive and 
correlated temporal behaviors. Regarding any weighted directed graph as being composed of the 
accumulation of paths between its nodes, our construction uses random walks of variable length to 
produce time-extended structures with adjustable features. The procedure is first described in a 
general framework. It is then illustrated in a case study inspired by a transportation system for 
which the resulting synthetic network is shown to accurately mimic the empirical phenomenology. 
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Many systems in nature or related to human activities 
are conveniently represented as networks of interacting 
units that can exchange material or information. This 
approach, combined with techniques from graph theory, 
statistical physics and data analysis, has led to countless 
interesting studies and insights in numerous fields [T]-[6]. 

Until recently, network structures were often regarded 
as being stationary, both for simplicity and because of 
scarceness of datasets on the time variability of connec- 
tivity patterns. Thanks to advanced acquisition tech- 
nologies and large scale production of time-resolved data, 
temporal information has become more accessible in nu- 
merous contexts, from communication networks [THTT] to 
proximity patterns [12^ il3j and infrastructure networks 
[14l|T5]. This has led to the recent surge of activity in the 
field of "temporal networks" [TF. Data analysis revealed 
the coexistence of statistically stationary properties and 
local variations, interaction burstiness [7H2Q] and the oc- 
currence of non-local repetitive patterns [151 [H] • These 
structural properties affect the dynamical processes tak- 
ing place on networks [22H3T] . Therefore, they have to 
be accounted for in modeling approaches. 

Despite the proliferation of temporally resolved 
datasets, strong limitations perdure. Indeed, data are 
often only accessible in restricted forms, such as single 
samples of limited sizes and statistical relevance. Com- 
parison of connection patterns at different times is not al- 
ways possible. Some datasets consist only of aggregated 
information and do not provide access to the temporal 
course of events. In such circumstances, it is necessary 
to be able to generate synthetic time-extended structures 
whose aggregation would reproduce the data at hand. 
This would enable one to go beyond the approximation 
of static networks and to incorporate dynamical compo- 
nents in network structures. 

Several models of temporal networks have been pro- 
posed in the literature l32ll36] . Their dynamics 
mostly consist of link updates and show that a global 




FIG. 1: (color online). Schematic representation of the con- 
struction procedure: a weighted directed graph Q can be re- 
garded as a superposition of paths. Unfolding these paths 
in time results in itineraries (of variable characteristics) that 
generate a temporal network A/". 



complex space-time organization can emerge as the re- 
sult of simple pairwise cooperation or competition rules 
between individual units. Dataset randomization is also 
employed to create null models against which the tempo- 
ral pattern complexity of empirical data can be evaluated 
[Ii[22j. 

The modeling approach developed in this Letter is to 
some extent more direct and pragmatic. It is intended to 
provide a versatile procedure for the construction of time- 
dependent networks with adjustable characteristics, that 
mimic bursty and possible repetitive behaviors, as well 
as extended temporal correlations, as they are observed 
in real world systems. The goal is to obtain realistic tem- 
poral structures independently of any access to datasets 
or assumptions about basic interaction mechanisms. 

The proposed procedure originates from the observa- 
tion that in many graphs representing systems such as 
transportation, trade or communication networks, con- 
nection patterns are the result of activity spreading 
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along noisy itineraries within a complex structure. The 
weighted directed graph used to represent this activ- 
ity is then given by the accumulation of such itinerary 
traces. Our construction proceeds in the opposite man- 
ner. Starting from a given weighted directed graph ^, it 
assumes that this graph is the superposition of random 
walks in a given time window W, and seeks to "unfold" 
these walks, as schematically illustrated in Fig. [l] 

Below, we first describe the construction in a general 
setting. We then apply it to a case study inspired by a 
transportation system, and use this real-world system to 
benchmark the synthetic network. Finally, we put for- 
ward a characterization of the topological and temporal 
correlations that takes into account the heterogeneity of 
the activity in the network. 

Temporal network construction. The construction 
takes as input a graph Q = (tS, S) and a time window 
W = {0, 1, • • • , T— 1}. Here, 5 is a set of nodes and ^ is a 
set of weighted directed edges, either obtained from data 
or built with prescribed properties such as in- and out- 
degree distributions and/or distributions of link weights 
(for definitions of degree, weight and strength, see e.g. 
[6|). The temporal network is built as follows. 

(i) First, the following random walk characteristics are 
chosen: distributions of starting times in W, of starting 
locations in 5, of walk lengths, and of residence times at 
each node. 

(ii) Then, the random walks are generated indepen- 
dently, one after another, using the characteristics spec- 
ified above. Each walk defines an itinerary that consists 
of a list of events {(io,ii,ti), (^1,^2,^2), • • • , {ii-i.H.U)} 
with ti < • • • < ti (Fig. [T|. At each step, the node ip is 
uniformly chosen among the out-neighbors of ip-i, and 
the residence time t^+i — tp is drawn from the residence 
time distribution associated with ip. Any walk reaching 
a node with no out-link terminates there, even if it has 
not yet reached its prescribed length. 

(iii) As the construction proceeds, the graph Q is con- 
tinuously altered as follows: if a walk passes through the 
link i ^ the weight Wij is decreased by 1; whenever 
Wij reaches 0, the edge i ^ j is discarded from the graph 
and cannot be used anymore. 

(iv) Walks are generated until all edges have been dis- 
carded from Q. Then, the process terminates. 

The temporal network JV is defined as the union of 
all the constructed itineraries. By construction, the 
set of weighted directed edges resulting from the collec- 
tion of all events in JV coincides with the original set 
£. Itineraries are interpreted differently depending upon 
context. For instance, in transportation systems, the 
event (i^j^t) is regarded as a material displacement from 
i to j at time t. In communication networks, events cor- 
respond to transmission of information. 

Note also that the construction can be used to model 
routine/seasonal processes occurring in consecutive time 
windows, using a "noisy deterministic" rule: for each n, a 
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FIG. 2: (color online). Prescribed (solid line) and realized 
(circles) distributions of random walk lengths in the tempo- 
ral network J\f. The discrepancy is due to the decay of the 
number of remaining edges in Q as the construction proceeds. 



tunable fraction of walks is redrawn in Wn = {tiT^ nT + 
1, • • • , (n + 1)T — 1}, while the remaining itineraries are 
repeated identically from Wn-i to Wn- 

We claim that, under suitable choices of input param- 
eters, the network M can mimic real-world features, in 
particular bursty temporal patterns and extended corre- 
lations. To support these assertions, we benchmark the 
synthesized network in a case study inspired by a cattle 
trade system analysis p!5l; note this example is aimed 
to illustrate the construction algorithm, rather than to 
accurately fit specific datasets. 

Case study. In this application, we consider \S\ — 10^ 
nodes representing farms. The graph Q is constructed 
using a variant of the uncorrelated configuration model 
[37. Following statistics reported in [15 , we use as in- 
and out-degree distributions the power-law P{kin,out) oc 
^in out (with cut-off at ^/\S\ to avoid degree correla- 
tions [37] )• Weights are also power-law distributed ac- 
cording to P{w) oc w~'^-^. Interpreting events in JV as 
material displacements, we impose flux conservation at 
every node so that in- and out-strengths balance, i.e. 
XI j ^ji — ^ij holds for almost all i G S (in prac- 
tice, more than 99% of S. 

We consider a temporal window W of length T = 10^ 
units ("days"), with periodic boundary conditions, to 
generate A/" from Q. To this aim, we use random walks 
with uniformly distributed starting times and locations. 
We assign each node i a residence time r^, drawn from 
P(r^) (X for 1 < < 60; a walk visiting node i stays 
for time r^, plus an additional random delay drawn from 
a Poisson distribution with mean (The power-law 

P(r^) is inspired by [TF; similar results are obtained for 
an exponential P{ri).) Walk lengths are generated using 
a Poisson distribution with average 10. Notice that since 
the number of graph edges decreases as the construction 
proceeds, long walks can seldom be achieved when few 
edges remain and the realized distribution accordingly 
exhibits an excess of short paths, see Fig. [2] 

Three main types of analysis have been proposed in or- 
der to study structural and temporal properties of a time- 
resolved transportation network [Ml [151 [21] : (i) statistics 
of aggregated networks on various time scales, (ii) distri- 
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butions of activity and inactivity periods and (iii) repeti- 
tion of connection patterns such as temporal paths. Here, 
we show that when apphed to the synthetic network A/", 
these diagnostics proffer characteristics very similar to 
those in empirical data. 
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FIG. 3: (color online). Distributions of node in-degrees, link 
weights, and node strengths for networks aggregated over in- 
tervals of various lengths. For integration lengths < 10^, the 
figures display distributions obtained by averaging over all 
corresponding intervals in W. Power laws (thick black lines) 
with expected exponents are plotted for reference. 



Real world data exhibit robust heterogeneous behav- 
iors at all temporal scales, as asserted by power-law 
distributions of in-degrees, strengths and link weights 
[121 [Ml [15]. These distributions remain almost stationary 
when integrated over distinct windows of equal length. 
As shown in Fig. |3j the same properties are observed in 
A/". These plots also validate the construction algorithm: 
statistics aggregated over the whole W reveal power-law 
distributions with slopes in close agreement with the cor- 
responding ones in Q. 




(a) Node activity time distrib. (b) Node inactivity time distrib. 



FIG. 4: (color online). Distributions of activity (a) and inac- 
tivity (b) periods aggregated over intervals of various lengths 
Tagg = 1,3, 10, 100, i.e. fractions of nodes that are active at 
least once in each of (resp. inactive during) n consecutive in- 
tervals of duration Tagg (the maximal value of n is T/Tg,gg). 
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FIG. 5: (color online). Left. Autocorrelation function C{t) 
vs. t, averaged over 40 realizations of A/", for two distinct 
walk length averages and compared to a null model where all 
links have been temporally randomized, equivalent to impos- 
ing that all walks in the construction have length 1. Right. 
Number of r-tolerant temporal paths (r = 5) occurring at 
least twice vs. their length, both in A/" and in the null model. 



Robust heterogeneous features in transportation net- 
works have also been observed in activity and inactivity 
period statistics, i.e. in the distributions of the periods 
during which a node (or a link) is continuously active or 
inactive. Such observables typically exhibit broad distri- 
butions [14 . Fig. |4] reports these distributions for nodes, 
obtained from A/" using various aggregation scales Tagg, 
from the most detailed ("daily", Tagg = 1) resolution to 
coarser timescales mimicking weekly or monthly aggre- 
gations (similar results are obtained for links). Activity 
distributions reveal power-law behaviors in quantitative 
agreement with empirical data [15 , with decay exponents 
in the same range of values (increasing from —6 to —3 as 
the aggregation interval Tagg increases). Moreover, distri- 
butions of inactivity periods are broader, more concave, 
and extend across all timescales, as in [15 . 

In addition to burstiness, the activity in A/" shows posi- 
tive temporal correlations, as can be anticipated from the 
nature of the construction. In particular, global activity 
correlations can be appreciated by using 



C{t) 
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where A{t) denotes the total number of events (i^j^t) on 
day t and (•)ti means time averaging over W. As Fig. [s] 
(left) shows, C{t) is positive for all t in an extended in- 
terval, whereas its analog in a null model obtained by 
reshuffling all links is for any t 0. Note that such a 
null model is equivalent to a temporal network construc- 
tion in which we impose all the random walks to have 
length 1. 

Beside global correlations, local structures can be high- 
lighted by exhibiting chronological link sequences that 
occur more often than by chance [151 EI] • A typical ex- 
ample is given by r-tolerant temporal paths, i.e. paths 
composed of consecutive active links within r days of 
each other. The number of such paths is much larger in 
empirical temporal networks than in reshuffled data [15] . 
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(Note that only the case r = 1 was considered in [E].) 
Fig. [5] (right), which shows the number of r-tolerant tem- 
poral paths that occur at least twice in A/", possesses sim- 
ilar features, providing further evidence of the presence 
of strong correlations in the synthetic network. 

Detailed activity correlations The nature of the tem- 
poral network construction, which relies on random 
walks, suggests that non-trivial local activity correlations 
should also simultaneously emerge at both topological 
and temporal levels. Indeed, if a node is active at a cer- 
tain time, then an elevation of activity should be detected 
in its topological and temporal neighborhood. (By con- 
trast, such correlations are nonexistent in a model where 
all links have been reshuffled.) Here, we provide quanti- 
tative estimates of these correlations. 

In heterogeneous networks such as the current syn- 
thetic A/", local activity patterns are extremely diverse 
and are hardly captured by global averages such as C{t). 
Therefore, prior to averaging, we first need to divide the 
nodes of S into categories according to their total activity 
(defined as the total number of events in W, starting from 
the considered node). The activity at a node in S ranges 
from 1 to over 200, with average |A/'|/|5| ~ 6. Accord- 
ingly, we define a node as "busy" if its integrated activity 
exceeds a given threshold, say e.g. 40 (results are qualita- 
tively insensitive to this choice). Furthermore, non-busy 
nodes are categorized into level sets of the distance d to 
the set of busy nodes S, where d is measured following 
the edges in S. Altogether, the categories S = = 0}, 
{d = 1}, . . . , {d = 4} account for more than 99% of S. 

For each node category C, the correlations between the 
activity in nodes of C at time t and in their neighborhood 
at time t -\- r are then measured by 

Fc,r{r) = (aB(i,r)(^ + ^)){iGC,tGW|sout(i,^)>0} 

where aE{t) = ^ Sout{j^t)/\E\ is the average activity 

jeE 

per node in a subset of nodes E at time t (so^tij^t) is 
the number of events starting from j at time t), B{i^r) 
is the ball of center i and radius r, and the average (•) is 
taken over alH G C and t eW such that 5out(^7 1) > 0. 

Fig. |6] displays "activity spike curves" Fc^r{T) vs. r 
for two categories, B and {d = 2}, and r = 1,2. In a 
null model where all links are reshuffled (equivalent to a 
model built with random walks of length 1), each such 
curve would be perfectly flat except at r = 0, where 
it would take a larger value. In contrast, each curve 
has a distinctive shape here that is almost insensitive to 
the network realization. Hence, for each category, these 
curves provide a characterization of activity correlation 
patterns in a topological and temporal neighborhood. 

The peaked shape of Fc,^(r) illustrates how the occur- 
rence of activity in the center of a ball is accompanied 
by an elevated activity in the ball for a number of days 
before and after. The peaks are substantially sharper 
for smaller balls (of radius 1), and are instead barely 
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FIG. 6: (color online). Fc,r{r) vs. r for 2 categories of nodes 
C = B (busy nodes) and {d = 2}, and two neighboring ball 
radii r = 1 and r = 2. Each panel displays the curves cor- 
responding to 10 realizations of A/". Dashed lines refer to the 
average activity 77 = \J\f\/T/\S\ per day per node. 



detectable when averaging in the ball of radius 2 around 
B. At large values of r, Fc^^(r) approaches a well-defined 
baseline level of activity. For r = 1, this baseline is higher 
around busy nodes than around = 2}, pointing to a 
higher level of mean activity in the neighborhood of S, 
as it can be anticipated. 

Concluding remarks. In this Letter, a simple frame- 
work for the construction of temporal networks display- 
ing bursty, repetitive and correlated behaviors has been 
presented, based on random walks on a predefined aggre- 
gated graph. The construction is sufficiently versatile to 
adapt to a wide range of situations, by appropriate tuning 
of input parameters. In particular, it would be interesting 
to apply it to initial graphs Q with community structure. 
As random walks tend to be trapped inside communities, 
the construction is expected to naturally give rise to the 
emergence of temporal communities. Another possible 
extension would be to use, instead of random walks as 
building blocks, spreading trees of propagation processes 
on networks. 

Our procedure enables one to obtain plausible and re- 
alistic instances of temporal networks when only the ag- 
gregated structure is known. Hence, it can be employed 
to compensate for the lack of time-resolved data or to 
provide alternative scenarios when access to empirical 
information is limited. Moreover, that the synthetic net- 
work has tunable properties is essential to assessing the 
influence of time-dependent features on the dynamical 
processes on networks, especially when only aggregated 
information is available. Finally, structural flexibility 
also makes it possible to test the relevance of various 
definitions of centrality measures for nodes and links in 
temporal networks. 
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